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ABSTRACT 

The relative motion of stars and other celestial objects in very wide pairs, separated 
by distances of the order of 1 pc, is strongly influenced by the tidal gravitational 
potential of the Galaxy. The Coriolis component of the horizontal tidal force in the 
rotating reference frame tends to disrupt such marginally bound pairs. However, even 
extremely wide pairs of bodies can be bound over intervals of time comparable to the 
Hubble time, under appropriate initial conditions. Here we show that for arbitrary 
chosen initial coordinates of a pair of stars, there exists a volume of the space of 
initial velocity components where the orbits remain bound in the planar tidal field 
for longer than 10 Gyr, even though the initial separation is well outside the Jacobi 
radius. The boundary of this phase space of stable orbits is fractal, and the motion 
at the boundary conditions is clearly chaotic. We found that the pairs may remain 
confined for several Gyr, and then suddenly disintegrate due to a particularly close 
rendezvous. By reversing such long-term stable orbits, we find that entrapment of 
unrelated stars into wide pairs is possible, but should be quite rare. Careful analysis 
of precision astrometry surveys revealed that extremely wide pairs of stars are present 
in significant numbers in the Galaxy. These results are expected to help discriminating 
the cases of genuine binarity and chance entrapment, and to make inroads in testing 
the limits of Newtonian gravitation. 
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1 INTRODUCTION 

The Galactic tidal force becomes comparable to the force 
of gravitational attraction of widely separated double stars 
at distances close to the Jacobi radius, which approximately 
equals 1.8 parsecs for twin stars of solar mass. Beyond the 
Jacobi radius, the planar component of the tidal force may 
cause a relatively rapid disruption of a wide pair. Stars 
emerging from disrupted binary systems or stellar associ- 
ations are stretched out by the Galactic potential into nar- 
row streams in the direction of Galactic rotation, a pro- 
cess, which is well de scribed by the epicycle approximation 
jMakarov et al.l [20041 ). This does not mean that all pairs of 
stars separated by distances greater than the Jacobi radius 
will be disrupted by the tidal force. The Coriolis component 
of force acting on a pair of stars in the co-rotating reference 
frame centered on the primary (i.e., larger mass) companion 
is proportional to the instantaneous relative ve locity of the 
secondary companion (|jiang fc Tremain3l2010l ). Under ap- 
propriate initial conditions, the Coriolis force can confine the 
relative motion of two stars to a finite volume of the phase 
space, even if their mutual distance is far greater than the 
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Jacobi radius. The resulting trajectories are similar to reg- 
ular orbits, in that they have a periodic character and the 
separation never exceeds a certain limit within the Hubble 
time. 

In this paper, we present and discuss the results of a 
large number of numerical integration of pairs of stars in 
the rotating, noninertial reference frame, performed with the 
J. Chambers' Mercury code, appropriately modified for this 
case (see Appendix). The aim of these numerical simulations 
is to demonstrate that a certain area of the phase space 
of initial parameters exists well outside the Jacobi radius, 
where the resulting orbits are regular, bound and stable. 
By retracing long-term stable orbits at the boundary of the 
stable area, we prove that accidental capture of two passing 
stars into a stable double system is possible in principle. 



2 STABILITY OF EXTREMELY WIDE 
DOUBLE STARS 

Fig. [1^ shows an example of a stable orbit for a hypothet- 
ical pair of stars with a total mass of 1.72 solar masses, 
traveling on a circular orbit at the orbital radius of the Sun 
around the Galaxy. The primary star is at the origin of the 
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reference frame (x = 0, y = 0), and the initial phase space 
parameters of the companion are Xq = 0, j/o = —4 • 10^ AU, 
xo = 4.8 • 10~^ AU d~^, yo = 0. The companion, moving in 
shifting double loops around the primary for 10 Gyr (which 
was our integration time) never wanders away by more than 
the initial separation. The Jacobi radius for this case is equal 
to 3.55- 10^ AU, assuming the values fo r the Ports constants 
deter mined from the Hipparcos data (|Makarov fc Murphvl 
l2007t ). Throughout this paper, the x axis is directed toward 
the Galactic center, the y axis toward the direction of Galac- 
tic rotation, and our computations are restricted to the pla- 
nar case, as explained in Appendix. 

Such long-term stable orbits at extreme separations are 
only possible when the orbital angular momentum is oppo- 
site to the Galactic rotation, i.e., when the relative motion 
is retrograde. As seen from the North Galactic pole, the 
orbital motion is always counterclockwise in stable pairs. 
This important property has been deduced for the classi- 
cal restricted three-body problem (|Henon|[l970l l and, more 
recently, for regular orbits of stars around open cluster s in 
the Galactic tidal potential (|Fukushige fc Heggi j [ioool ) . A 
co-rotating pair of stars beyond the Jacobi radius is dis- 
rupted by the Coriolis acceleration well before it completes 
a full revolution in the Galaxy (approximately, 230 Myr). 
Thus, the chances of survival of very wide double systems 
are defined by initial conditions, barring external perturba- 
tions from other field stars. By integrating numerous tra- 
jectories at a step of 0.1 AU d^^ in xq and yo for the pre- 
viously used initial coordinates xq — 0, yo = —4 • 10^ AU, 
we found that the area of stable orbits is approximately 
triangular in shape defined by the vertices at {xo,yo) = 
{(-4.74, 0), (-h4.74, 0), (0, 7.0)} • 10"^ AU d'^ Most of the 
trajectories with initial velocities inside this triangle are 
bound in the long term (over 10 Gyr or longer). Most of 
the trajectories outside this triangle are not bound, with the 
distance rapidly increasing beyond the initial separation. 



Figure 2. Survival time of pairs of stars initially separated by 
400 000 AU computed with a step of 0.17 mm s~^ in initial ve- 
locity components. A tiny variation in phase space parameters 
may cause a drastic change in stability at this transitional region 
separating bound and unbound trajectories. 



3 CHAOS AND ENTRAPMENT 

The boundary of the phase space of stable orbits assisted 
by the Galactic tidal field is fractal, as established numeri- 
cally by performing integrations at a small step in the initial 
conditions. In the vicinity of this boundary, a tiny change 
in the initial velocity, for example, may alter the fate of the 
stellar pair, resulting in orbits that are bound for longer 
than 10 Gyr, or ones that fall apart within 100 Myr. Fig. 
[2] shows the survival time of pairs with the initial coordi- 
nates at Xo = 0, yo — —4 • 10^ AU, and initial velocities 
varied by small amounts around io — 4.7394- 10~^ AU d~^, 
yo = 0. Sudden changes in survival time take place on scales 
as small as 0.17 mm s~^ in either dimension of the velocity 
space. In a separate experiment, we obtained a stable orbit 
over 10 Gyr at io = 4.7398 • 10"^ AU d"\ an orbit which 
disintegrated in 3.26 Gyr at io = 4.73981 • 10"^ AU d"\ 
and, most remarkably, an orbit which fell apart after 9 Gyr 
at io = 4.7398000004-10"^ AU d"^ The machine precision 
of these computations is approximately one part in 10^^ . 

The motion of stars with initial parameters at the 
boundary of the space of stable orbits is strongly chaotic. 
Using the method to estima te the Lyapunov time suggested 
for the outer Solar System (|Haves 2007, 2008), two sibling 
trajectories were integrated at xo = 0, yo = —4 • 10^ AU, 
io = 4.7394 ■ lO"'^ AU d"\ yo = by changing the initial 
value io by one part in 10^^. The spatial separation between 
the siblings grew exponentially over the first few Gyr, with 
an overlaid oscillation of 0.72 Gyr period. The estimated 
Lyapunov time was 0.56 Gyr, which is approximately equal 
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Figure 3. Trajectory of a 0.72 Mq companion around a 1 Mq 
star in the plane of the Galaxy, entrapped into a bound pair for 
8.3 Gyr by the Galactic tidal force. Both the capture and the 
ejection of the companion were achieved through gravitational 
maneuvers in the vicinity of the primary component. 



to the duration of two Galactic revolutions. Both orbits were 
found stable over 10 Gyr. 

The long-term stable obit at xo = 4.7398000004 ■ 10"^ 
AU d~^ , discussed above, is of special interest for the prob- 
lem of capture of random passing stars into double systems. 
This pair was bound for 9 Gyr, after which it suddenly fell 
apart when the secondary component wandered too close to 
the primary and was accelerated by its gravitational pull. 
Once the pair is torn apart, the distance between the com- 
ponents begins to grow rapidly. We can construct a long- 
term stable pair by taking any of the points on the outgoing 
trajectory, reversing the time and performing a reflection 
in the four- dimensional phase space. The reflection is neces- 
sary, because only retrograde trajectories are stable. A stable 
trajectory obtained this way for a pair of initially unbound 
stars is shown in Fig. O A chance fly-by star was decel- 
erated by the tidal force, went through a close periastron 
where it lost more kinetic energy and started to orbit the 
primary on the characteristic double loops we have seen in 
Fig. [T] After ~8 Gyr of quasi-regular motion, the companion 
wandered too close to the primary and was ejected again. 
Obviously, the original trajectory, which was bound for ~9 
Gyr, could not be exactly reproduced in this experiment. 
The chaotic motion and the fractal character of the bound- 
ary between long-term and short-term trajectories make the 
machine precision (about 16 significant digits) of our cal- 
culations insufficient to match two orbits over such a long 
time. Thus, accidental entrapment of passing stars into sta- 
ble pairs is possible in principal, but should occur extremely 
seldom because of the small area of suitable phase space pa- 
rameters. Furthermore, any perturbation from a third pass- 
ing body will drastically change the fate of the system, in 
most cases resulting in a relatively rapid breakup. Occa- 
sionally, the disturbance from a chance fly-by will drive an 
initially unstable pair into the space of stable trajectories, 
but the probability of such events at separations above the 
Jacobi radius is very small. This entrapment mechanism is 
similar to the chaos-assisted capture of irregular moons by 
the Sun- Jupiter system, which can only happen in the thin 



layer of chaot ic motion s eparating the regions of scattering 
and stability (jAstakhov et al.ii2003 '). 



4 DISCUSSION 

It follows from these considerations that the entrapment of 
random stars into extremely wide systems has low probabil- 
ity and that such associations should be short-lived. Yet, a 
growing body of observational evidence indicates that pairs 
of stars with separations up to 0.1 - 1 pc are abundant even 
in the close neighborhood of the Sun l|Makarov et al. I l2008l : 
ICaballerdl2009l : IShava fc Olline"201lV Some of these pairs 
may be the r emnants of disintegrating open clusters or OB 
associations (|Kouwenhoven et al1l2010l ). but not all of the 
components are sufficiently young. If verified by accurate 
astrometric observations, the existence of only a few pairs 
of stars with separations greater than the Jacobi radius in 
the Solar neighborhood will challenge the classical theory 
of gravitation in the domain of weak forces. Any modifica- 
tion of Newtonian dynamics in the weak-field regime that 
boosts the gravitational acceleration at large distances will 
drastically increase both the probability of entrapment of 
random pairs of field stars and the survival time of very 
wide binaries. Even the nearest known stellar system to 
us, the Alpha Centauri AB and the remote Proxima Cen- 
tauri companion, may pose a problem for Newtonian dy- 
namics because of the relatively large difference in v elocity 
/An osova et al. 1 1 19941 : IWertheimer fc Laughlinll2006l ). iBeechI 
(200^) suggested that the apparent stability of t his triple 
syste m may provide a crucial test for the MOND (jMilgroml 
Il983i ). or the theory of a coupling betwee n the scalar cur- 
vatur e and the matter Lagrangian density (|Bertolami et al.l 
|2003), where an additional force can widen the allowable 
limit on the orbital acceleration. For the case studied in this 
paper (Mtot = 1.72Mq, j/q = -4 • 10^ AU), the Newtonian 
gravitational acceleration is much smaller than the thresh- 
old acceleration of MOND (ao ~ 1.2 ■ 10""^ m s~^). Using 
the asymptotic expression for acceleration of the companion 
star, a — y/G Mtot ao , we estimate the distance of the 
Lagrange points to be 

/ VGA/tot ao\ ' 

4An ) 

which for this pair of stars is rj = 2.2 • 10^ AU, or 10.5 pc. 
The Jacobi radius in MOND is 6.2 times larger than the 
Newtonian Jacobi radius, and the area of stable orbits is 
roughly 38 times greater in the {x, y) plane. Very accurate 
radial velocity measuremen ts of Proxim a Cen would be re- 
quired for a conclusive test ([Beec^lloiJ). Very wide pairs of 
nearby stars, due to even stricter criteria of survival, may 
provide the most sensitive experimental data for testing the 
theories of gravitation in the weak-field regime. 
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APPENDIX A: GALACTIC TIDAL FORCE 

The planar components of the Galactic tidal accelera- 
tion are accurately approximated by the epicycle model 

for all objects moving on nearly c i rcular orbi ts at 

the solar radiusllJiang & Tremaind l2010l : iKing et"aL 1 1 19901 : 
iMakarov et al.ll2004^ : 

X = 4:Anx - 2Q.y 
y = 2nx 

where A is the Oorts constant, and fl is the rotation fre- 
quency. The vertical component of acceleration, to a good 
approximati on for small vertica l velocities at the plane, 
is harmonic (|Makarov et al.ll2004l ). and, therefore, tends to 
compress very wide binaries. It is sufficient to consider 
the stability problem restricted to the planar case, ignor- 
ing the vertical dimension. All inte grations were performed 
with the well-tested Mercury code (|Chambers & Migliorinil 
1 19971 ) by adding the tidal acceleration components in the 
user-defined external acceleration subroutine. In the units 
adopted in the Mercury code, the assumed parameters were 
AAQ. = 1.19 • 10"^° d-2 and 2^ = 1.5 • 10"^° d"^ We used 
the Bulirsch-Stoer option of integration and set the integra- 
tion step to 200 d. The companion was considered ejected 
when the distance from the primary exceeded 900 000 AU. 
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